Non-monotonic behavior of multiplicity fluctuations 
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Abstract. We discuss recently measured event-by-event multiplicity fluctuations in relativistic 
heavy-ion collisions. It is shown that the non-monotonic behavior of the multiplicity fluctuations 
as a function of collision centrality can be fully explained by the correlations between produced 
particles. 



1. Introduction 



> 

cn 
(N 
O 

00 ' Event-by-event fluctuations in heavy-ion collisions have been recently measured both at CERN 

. SPS and BNL RHIC (see a brief review ^). The data (usually used F or ^ measures to quantify 

the fluctuations of transverse momentum) which show a non-trivial behavior as a function of 
collision centrality, have been theoretically discussed from very different points of view, including 
complete or partial equilibration, critical phenomena, string or cluster percolation, production 
of jets. In spite of these efforts a mechanism responsible for the fluctuations is far from being 
^ ■ uniquely identified. 

^1 In [2] we have established connection of the ^-measure with the fluctuations of multiplicity 

and advocated that fluctuations of transverse momentum are practically irrelevant for the 
considered problem, i.e., for the apparent structure seen in collision centrality dependence. 
rS I We argue [3] that the non-trivial behavior of the so called transverse momentum fluctuations 

■ (quantified by F or <I> - measure) can be fully explained by the multiplicity fluctuations. 

Recently, the NA49 Collaboration published very first data on multiplicity fluctuations as 
a function of collisions centrality jlj, see also in this volume. Unexpectedly, the ratio 
Var{N)/ {N) , where Var{N) is the variance and (A^) is the average multiplicity of negative 
particles, changes non-monotonically when number of wounded nucleons grows. It is close 
to unity at fully peripheral {Ny^ < 10) and completely central {Nu, > 250) collisions but it 
manifests a prominent peak at ~ 70. The measurement has been performed at the collision 
energy 158 AGeV in the transverse momentum and pion rapidity intervals (0.005, 1.5) GeV and 
(4.0,5.5), respectively. The azimuthal acceptance has been also limited, and about 20% of all 
produced negative particles have been used in the analysis. 

The aim of this paper is to show that the observed effect (non-monotonic dependence of the 
ratio Var{N)/ < N > on the number of projectile participants Np, i.e., on the production 
volume V) stem from the correlations between produced particles and the correlations are of 
nuclear and electromagnetic interactions origin. 



2. Correlations and fluctuations 

Let us introduce the single-particle configuration distribution function jSj 

ni(ri)=n (1) 
where n is (constant) density of particles. Two-particle configuration distribution function: 

n2(ri,r2) = • n2(|ri - r2|) = • n2(r) (2) 
is expressed by two-particle correlation function j^2(?') following: 

n2{r) = I + U2{r) (3) 

where, asymptotically for r oo, i'2if) = and n2{r) = 1. 
For the volume V one gets 

< A^y >= J (irini(ri) = J drin (4) 
V V 

<Nl>= j (irini(ri) + J dri J (ir2n2(ri, rs) (5) 



V V 



The variance equals: 

< N^> - <Nv = J drini{ri) 



+ j dri y"dr2[n2(ri,r2) - ni(ri)ni(r2)] (6) 



V V 



and the normalized variance is given as: 



< Af2 > _ < >2 

— = 1 + / or [n • n2(rj — nj (7j 



< A^y > 

Substituting now Eq. (jSJ to Eq. Q one imediately leads to the formula: 

< iV^ > _ < ^ >2 



< iV > 



l + n dru2{r) (8) 



which shows clearly that the normalized variance is given by two-particle correlation function. 
Because: 

VariN) 

—^ = l + nV <V2> (9) 

where we note that the integral / dri / dr2i^2(|ri — r2|) =< V2 > V^, then it is clear that: 

V V 

( VariN) _ ,\ 

< -2 >= (10) 

has numerical values close to zero. Determined this way < f2 > do not depends on the detector 
acceptance p. Because, the normalized variance for the accepted particles: 



Var{N) 

<N > 



[l-p)+p 



(11) 



and the accepted multiplicity: 



<N>=p-<Np=i> (12) 

one finds from Eq. (|1U() the correlation function for the accepted particles < 1^2 > = < 1^2 >p=i- 
Fig. n shows mean value of correlation function < 1/2 > calculated from data on the ratio 
Var (N) / < N > obtained at 158 AG eV ^. 
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Figure 1. Mean value of correlation function 
calculated from data obtained at 158 AGeV . 



Figure 2. The dependence ol P ■ W (r) on 
relative distance r. 



2.1. Simple example 

Let the correlation function be in the form: 



In this case: 



V2{vi - r2) = V2{r) = exp(-2r/A) 



(13) 



Var (N) 
<N> 



1+ < N > / / dri(ir2i/2(ri - r2) 



V V 
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1+ < > {R- r)u2{r)dr 



(14) 



and one gets 



<N > 

It the limit R >> X we have 



Var(N) <A^>/A\V / R\ R 



(15) 



(16) 



then for R « \ 

Var{N) 
<N > 



1+<N> (17) 



3. Correlations and interactions 

Two-particle distribution function n2{r) normalized for mean density, in the equilibrium case 
(with inverse temperature /3) is given by the following Boltzmann factor [H]: 

n-^n2{r) = n-exp{-(3Wir)) (18) 
For the rigid spheres with the radius dQ, where 

mr)= T '^o? (19) 
^ ^ r > 2(^0 



we have the correlation function: 

n2{r) = r < 2do 
V2{r) = -1 r < 2do 

3.1. Two potentials 

Let us consider two potentials: the electrostatic Debye potential [Zj 



(20) 



exp j-r/Xe) 

(fe{r) = e (21) 

r 

and the nuclear Yukawa potential ISj: 

exp j-r/Xn) 

fn{r) = -g (22) 

r 

The attractive (nuclear) and repulsive (electrostatic) interactions leads to the effective energy: 

W{r) = e-y^e{r)-g- (r) (23) 

The dependence of f3-W (r) on relative distance r is presented in Fig. [21for the numerical values: 
Xn = 1.5 /m] and Ae = 3.5 fm. 

For the radius r > 2do one can write the correlation function 

/X / \ ^ ( exp(-r/A„) exp(-r/Ae)\ , . 

V2[r) = n2(r) - 1 = exp a„ ae - 1 (24) 

\ r r J 

where strength factors are equal to: an = Pg"^, ae = /?e^ = (47rnAg) ^ and interaction lengths 

are: A„ = h{2TTmc)~^, Ag = (47re^n/3) 
In the first approximation one can write 

f \ ^^P i-^/^n) exp (-r/Ae) . . 

(r) PS a„ ae (25) 

r r 

what allows for analytical integration n f (irz^2(?') and analytical estimation of normalized 
variance Var{N)/ < N >. 

The problem of local and global fluctuations (the so called finite-size effect) was originally 
discussed by Nicolas [0] for correlation function: 



.2(r)=a^^P(^ (26) 
r 

Integrating it (via the Fourier transformations method) in the spherical volume V with the 
radius R gives [TO] : 



r^r. ^ ^ Var(N) ^ f , exp(-r/A) 
f{R,X,a) = =l + nj dra ' ' 



A5 

1 + 67ran— -7 
R-^ 



A +3 A -exp(-2fl/A)^ 1 + - 



(27) 



Asymptotically for small R {V = A/'inR^ « A'^) one gets: 



/ (i?, A, a) = — — — ocl + (28) 

< N > 5 

whereas for large R {V = 4/37ri?^ >> A'^) one gets: 

/WA,a) = ^o<W(l-fA) ,29) 

For our correlation function 1^2 (1^), given by Eq. (|25|) . one immediately gets: 

Var (N) 



<N > 



f{R,Xn,an)-fiR,Xe,ae) (30) 



In this paper the integral Eq. (jH)) has been calculated via Monte-Carlo method, assuming 
the homogenous spherical production volume, V ~ Np, to be proportional to the number of 
participants Np. Fig 13] shows normalized variance of the multiplicity distribution and mean 
multiplicity of negatively charged particles produced in collisions at 158 AGeV fitted by the 
two potentials model. It is interesting to note that different shapes of the production region 
(with particle density distribution n (r)) influence the presented results for the large number of 
participants Np and allow us to obtain much more satisfactory fit to experimental data. 

3.2. Single potential (dipole-dipole interactions) 

Electrostatic correlations play an important role in analysis of thermodynamic system jlll I12j . 
In following, we analyze one more (speculative, but possible) scenario, which result in non- 
monotonic behavior of fluctuations. 

Let us assume that the particles are produced in pairs (vr+vr^) creating dipoles. The potential 
describing interaction between pi and p2 dipoles is given by: 

. ^ ^ ( 3(pir) (p2r) ^ 
'^W = ^[P^P^ -2 J 

< p >^ 

= g — (sin 6i sin 62 cos (</>i — ^2) + cos 9i cos ^2 — 3 cos 61 cos 62) 

= ^73^/ (^''^) (31) 



and changes from —2 < p >^ /r^ to +2 < p /r^, depending on dipoles space arrangement. 
Changing the polarization of dipoles depending on distance between them / {9, (j)) ^ f (r, 9, 4>) 



allows for transition from negative values of potential (low values of r; particles correlation) to 
positive (big r; particles anti-correlation). For instance the dependence of the type: 



/(r,0,</.) = -2exp 



{r/xy 



+ 2 1 - exp 



{r/xy 



(32) 



provides to the correlation function i^2('") which allow for approximate description of 
experimental data (c.f. Fig^J where we fit data with the parameter A = 2.6 /m). One can 
get the effect of changing / (r, 9, (j)) via correlation between the dipole angles, for example: 



= TT ( 1 — exp 



ir/Xf 



(33) 




Figure 3. Normalized variance of the multi- 
plicity distribution (a) and mean multiplicity 
(b) of negatively charged particles produced 
in collisions at IbSAGeV |4j fitted by the two 
potentials model. 



Figure 4. Normalized variance of the multi- 
plicity distribution (a) and mean multiplicity 
(b) of negatively charged particles produced 
in collisions at Ih^AGeV [Ij fitted by the 
dipole model. 



4. Connections to percolation 

The particles situated in correlation length A distance each other builds the cluster. For low 
R < A/2 the whole production volume belongs to the single cluster. If the size of the fireball 
increases "the percolation" to the next clusters appear jl^j. The number of percolation clusters 
Mc is proportional to: 



(34) 



The number of particle pairs in distance d, included in the sphere with radius R is distributed 
as: 



d 
R 

2 - 



A 
R 



0<d< R 
R<d<2R 



(35) 



and the number of particles pairs with relative distance d < A to the total number of pairs is 
given by: 



2R 



Mp 



1 



R> X 

^<R<X 

R<1 



(36) 



In the simplest approximation the number of particle pairs for single cluster equals Mp /Mq and 



correlation function V2 ~ {Mp/Mc)^ is given by: 



R>\ 



^ <R<X 



R<^ 



(37) 



5. Connections to nonextensivity 

Let us now continue discussion of the multiplicity fluctuations from the point of view of 
its possible connection with nonextensivity (for other hints on nonextensivity in hadronic 
production processes and references to nonextensive statistics, see [21 E])- When there are 
only statistical fluctuations in the hadronizing system one should expect the Poissonian form 
of the multiplicity distribution. The existence of intrinsic (dynamical) fluctuations result in 
broader distributions usually expressed via Negative Binomial form and characterized by the 
parameter (so called the inverse number of " clans" ) : 



k- 



Var{N) 



1 



<iV>2 <iV>' 
By using Eq. ©, Var (N) / < N >= 1+ < N >< U2 >, then one gets: 



q 



1 = k 



-1 



= < 1^2 > 



(38) 



(39) 



where q is an nonextensivity parameter. 

Small systems (with diameters non exceeding correlation length) are characterized by the 
parameter q > 1 and just in the thermodynamic limit q = 1 (for big size systems with 
R >> \, where the correlation effects are small and all subsystems are characterized by the 
same temperature). 

The q parameter could be connected with the number of correlation clusters ("copies" of the 
system) Mq- If Mc grows, the Tsallis entropy pTS] : 



S = lim 



1 



where nonextensivity parameter q 
entropy: 



g-i 1- 
1 + m ~ 1 + 



Mc 



(40) 

leads to the thermodynamic Shannon 



S = — y^plnp 



(41) 



because: 



Inp = lim — 

m— >0 

Experimental data prefers dependence of the type: 



(42) 



m 



Mc ~ exp 




) 



(43) 



6. Conclusions 

The observed effect (non-monotonic dependence of the ratio Var(N)/ < N > on production 
volume V ~ Np) stem from the correlations between produced particles. The correlations are 
of nuclear and electromagnetic interactions origin. It is possible to describe the appearance of a 
prominent peak in Var{N)/ < N > at number of participants Np ~ 35, for charged particles. 
For the neutral particles we observe lack of non-monotonic behavior of Var{N)/ < N > when 
number of wounded nucleons grows (satisfy the experimental observation J7j). Effect depends 
on energy, because particles density n changes with energy (the magnitude of the effect ~ n and 
position of maximum also shifts, because correlation length Ae — l/^/n). 
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